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Based  on  the  coupled  density  and  energy  balance  equations,  a  dynamical  model  is  proposed  for  exploring 
many-body  effects  on  optical  carrier  cooling  (not  lattice  cooling)  in  steady  state  in  comparison  with  the  earlier 
findings  of  current-driven  carrier  cooling  in  doped  semiconductors  [X.  L.  Lei  and  C.  S.  Ting,  Phys.  Rev.  B  32, 
1112  (1985)]  and  tunneling-driven  carrier  cooling  through  discrete  levels  of  a  quantum  dot  [H.  L.  Edwards  et 
al.,  Phys.  Rev.  B  52,  5714  (1995)].  This  dynamical  carrier-cooling  process  is  mediated  by  a  photoinduced 
nonthermal  electron-hole  composite  plasma  in  an  intrinsic  semiconductor  under  a  thermal  contact  with  a 
low-temperature  external  heat  bath,  which  is  a  generalization  of  the  previous  theory  for  a  thermal  electron-hole 
plasma  [H.  Haug  and  S.  Schmitt-Rink,  J.  Opt.  Soc.  Am.  B  2,  1135  (1985)].  The  important  roles  played  by  the 
many-body  effects  such  as  band-gap  renormalization,  screening,  and  excitonic  interaction  are  fully  included 
and  analyzed  by  calculating  the  optical-absorption  coefficient,  spontaneous  emission  spectrum,  and  thermal- 
energy  exchange  through  carrier-phonon  scattering.  Both  the  optical  carrier  cooling  and  heating  are  found  with 
increasing  pump-laser  intensity  when  the  laser  photon  energy  is  set  below  and  above  the  band  gap  of  an 
intrinsic  semiconductor.  In  addition,  the  switching  from  carrier  cooling  to  carrier  heating  is  predicted  when  the 
frequency  detuning  of  a  pump  laser  changes  from  below  the  band  gap  to  above  the  band  gap. 
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I.  INTRODUCTION 

It  is  well  known  that  the  optical  properties  of  an  atomic 
vapor  are  essentially  determined  by  the  properties  of  a  single 
atom  in  the  linear-response  regime,1  i.e.,  a  single-particle 
problem.  However,  the  photoexcited  electrons  and  holes  in¬ 
side  an  undoped  semiconductor  are  quasiparticles  moving 
freely  throughout  the  whole  crystal  in  the  absence  of  lattice 
vibrations  and  structural  defects.  Therefore,  the  long-range 
characteristics  of  the  Coulomb  interaction  between  charged 
particles  lead  to  a  genuine  many-body  problem  to  determine 
optical  properties  of  a  semiconductor.  The  strong  Coulomb 
interaction  between  electrons  or  between  holes  forces  these 
charged  quasiparticles  to  constitute  a  plasma,  while  the 
strong  Coulomb  interaction  between  an  electron  and  a  hole 
forces  the  pair  to  form  an  exciton.  The  many-body  effects 
cause  laser-intensity  (or  charge  density  which  is  proportional 
to  the  laser  intensity)  dependence  of  the  optical  spectra  for 
semiconductors  in  the  linear-response  regime,  i.e.,  the  band- 
edge  optical  nonlinearity. 2-4 

When  an  n-doped  semiconductor  is  under  a  bias  field  be¬ 
tween  two  electrodes,  free  electrons  in  the  semiconductor 
undergo  a  center-of-mass  motion  with  a  drift  velocity.5  In 
this  situation,  free  electrons  reach  a  quasiequilibrium  state 
with  a  resulting  electron  temperature  due  to  an  ultrafast  rela¬ 
tive  scattering  motion  of  electrons  with  themselves,  with  lat¬ 
tice  ions  and  impurities.  However,  this  electronic  system 
usually  does  not  stay  in  a  thermal  equilibrium  with  the  lattice 
that  has  a  different  constant  heat-bath  temperature.  There¬ 
fore,  it  is  possible  for  drifting  electrons  to  transfer  their  ther¬ 
mal  energy  to  the  lattice  even  when  they  are  at  a  lower  tem¬ 
perature  than  that  of  the  lattice  due  to  reduced  phonon  energy 
induced  by  a  Doppler  shift.5  As  a  result,  electrical  carrier 
cooling  in  steady  state  under  a  bias  field  is  expected. 

Another  proposed  approach  for  electrical  carrier  cooling 
involves  using  a  tunneling  band  structure  in  quantum  dots,6 


in  which  hot  electrons  above  the  electron  Fermi  energy  (or 
hot  holes  below  the  hole  Fermi  energy)  are  extracted  using 
selective  resonant  tunneling  through  a  discrete  energy  level 
of  a  quantum  dot.  This  tunneling  process  modifies  the  qua¬ 
siequilibrium  carrier  distribution  function,  leading  to  a  sharp¬ 
ened  step  feature  in  a  distribution  (a  downward  step  at  the 
Fermi  energy)  with  a  lower  temperature. 

On  the  other  hand,  when  an  undoped  semiconductor  is 
subjected  to  a  strong  optical  pumping,  photogenerated 
electron-hole  ( e-h )  plasmas  are  usually  not  in  thermal  equi¬ 
librium  with  a  lattice  having  a  fixed  heat-bath  temperature. 
This  issue  has  not  been  addressed  in  the  previous  research,2-4 
although  it  may  be  crucial  for  studying  the  saturation  of  op¬ 
tical  absorption  in  undoped  semiconductors  under  strong  op¬ 
tical  pumping.  When  an  absorption  saturation  is  reached  in  a 
steady  state,  there  still  exists  an  energy  imbalance  between 
the  optically  absorbed  power  and  the  power  loss  due  to  spon¬ 
taneous  emission  of  photons.  Therefore,  a  dynamical  energy 
equation  is  required  to  describe  this  energy  imbalance  and 
the  compensation  from  the  thermal  exchange  between 
phonons  and  photoexcited  carriers.7  As  a  result,  optical  car¬ 
rier  cooling  is  expected  in  steady  state  with  a  pump  laser  if 
the  optically  absorbed  power  becomes  smaller  than  the  spon¬ 
taneous  power  loss,  which  is  different  from  the  lattice  cool¬ 
ing  in  semiconductors  by  photoluminescence.  This  turns  out 
to  be  the  microscopic  origin  for  laser  cooling  of  a  lattice  in 
thermally  isolated  semiconductors7-10  if  the  heating  to  the 
lattice  such  as  Auger  recombination11  can  be  controlled  in 
the  system. 

Carrier  cooling  and  lattice  cooling  are  physically  differ¬ 
ent,  although  they  may  possess  some  mathematical  similari¬ 
ties.  The  latter  is  a  net  cooling  with  a  thermal  isolation  of  the 
whole  system  from  its  environment,  while  the  former  is  only 
a  partial  cooling  of  a  system  with  a  lower  carrier  temperature 
than  the  lattice  temperature  set  by  the  thermal  contact  to  an 
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external  heat  bath.  In  addition,  the  very  difficult  lattice  cool¬ 
ing  in  semiconductors,  if  it  exists,  could  be  observed  from 
the  shift  of  the  photoluminescence  peak,  while  the  easy  car¬ 
rier  cooling  can  be  directly  verified  by  the  delicate  tempera¬ 
ture  dependence  of  the  resonant  tunneling  current.  Concep¬ 
tually,  slow  lattice  cooling  with  quasiequilibrium  distri¬ 
butions  of  carriers  and  phonons,  as  well  as  an  adiabatic 
variation  of  the  lattice  temperature,  requires  different  carrier 
and  lattice  temperatures  so  that  the  thermal  exchange  be¬ 
tween  carriers  and  phonons  becomes  possible.  Therefore,  the 
basic  assumption  of  equal  temperatures  in  the  papers  by 
Rupper  et  al.12  and  by  Sheik-Bahae  et  al. 13  for  both  carriers 
and  phonons  cannot  be  applied  to  quasiequilibrium  distribu¬ 
tions  of  earners  and  phonons. 

In  this  paper,  we  will  introduce  the  coupled  dynamical 
density  and  energy  balance  equations  to  search  appropriate 
laser  photon  energies  and  intensities  for  demonstrating  the 
steady-state  optical  carrier  cooling  mediated  by  nonthermal 
electron-hole  plasmas  in  an  intrinsic  semiconductor  with  a 
fixed  lattice  temperature  set  by  an  external  heat  bath.  We  will 
include  many-body  effects  such  as  band-gap  renormaliza¬ 
tion,  screening,  and  excitonic  interaction  in  calculating  the 
optical-absorption  coefficient,  spontaneous-emission  spec¬ 
trum,  and  thermal-energy  exchange  through  carrier-phonon 
scattering.  By  adjusting  the  laser  photon  energy  above  (be¬ 
low)  the  band-gap  energy  of  an  intrinsic  semiconductor,  we 
look  for  optical  carrier  heating  (cooling)  in  the  system.  In 
addition,  by  sweeping  the  frequency  detuning  of  a  pump  la¬ 
ser  with  respect  to  the  band-gap  energy,  we  expect  to  see  a 
switching  from  optical  carrier  cooling  to  carrier  heating. 

The  outline  of  this  paper  is  as  follows.  In  Sec.  II,  by 
employing  the  coupled  dynamical  density  and  energy  bal¬ 
ance  equations,  we  present  our  model  and  theory  to  explore 
many-body  effects  on  optical  carrier  cooling  mediated  by  a 
nonthermal  e-h  composite  plasma  in  intrinsic  semiconduc¬ 
tors.  In  Sec.  Ill,  we  display  and  discuss  numerical  results  for 
the  optical  carrier  cooling  as  function  of  the  laser  intensity 
and  photon  energy,  demonstrating  the  important  role  played 
by  the  many-body  effect  in  the  process  of  optical  carrier 
cooling.  A  brief  conclusion  is  given  in  Sec.  IV  with  a  remark. 

II.  MODEL  AND  THEORY 

The  ultrafast  intraband  electron-electron  (hole-hole)  scat¬ 
tering  in  an  electron  (hole)  plasma  ensures  all  electrons 
(holes)  in  a  quasiequilibrium  state  with  a  time-dependent 
electron  (hole)  temperature  Te(t)  [Th(t)]  through  an  ultrafast 
thermalization  process.14  Furthermore,  the  strong  interband 
Coulomb  interaction  between  the  electron  and  hole  plasmas 
locks  their  temperatures  to  the  same  value  Tc=Te=Th  (Ref. 
15).  The  lattice  system  is  assumed  in  thermal  contact  with  an 
external  heat  bath  at  a  fixed  temperature  T L.  As  a  result,  the 
phonons  with  a  different  temperature  usually  do  not  reach  a 
thermal  equilibrium  with  the  charged  carriers  in  the  plasma.7 

A.  Dynamical  density  equation 

The  relative  importance  of  various  scattering  processes  in 
an  undoped  semiconductor  depends  on  the  type  of  materials 


and  the  quality  of  material  growth.  For  commonly  used  III-V 
materials  such  as  GaAs,  the  growth  processes  are  usually 
highly  controlled,  such  that  the  density  of  crystal  defects  is 
low  enough  to  become  negligible.  In  this  paper,  we  only 
consider  the  intrinsic  properties  of  semiconductors  for  pos¬ 
sible  optical  carrier  cooling,  i.e.,  Tc  <  7)  .  Therefore,  a  defect- 
free  crystal  is  assumed  for  the  host  semiconductor  to  be  con¬ 
sidered.  Based  on  the  conservation  of  number  of  electrons  Ne 
and  holes  Nh  per  unit  volume,  the  dynamical  density  equa¬ 
tion  that  describes  the  temporal  changes  of  the  density  Nc 
=Ne=Nh  of  photoexcited  carriers  in  intrinsic  (defect-free  and 
undoped)  semiconductors  under  a  spatially  homogeneous  la¬ 
ser  illumination  can  be  written  as3 

dNc  /3abs(flz.)/o  _  .  . 

dt  ~  mL  sp’ 

Here,  hil,  in  Eq.  (1)  is  the  energy  of  incident  photons,  70  is 
the  energy  flux  of  the  incident  laser  beam,  /3abs(07  )  is  the 
linear  optical-absorption  coefficient  of  the  semiconductor, 
and  lZsp  is  the  total  rate  of  spontaneous  emission  of  photons 
per  unit  volume.  The  three-body  Auger  recombination  of  ex¬ 
cited  e-h  pairs,11  which  is  proportional  to  the  cube  of  Nc,  has 
been  neglected  in  Eq.  (1)  for  low  densities  Nc  and  low  laser 
intensity  70.  The  steady-state  solution  of  Eq.  (1)  predicts  a 
saturation  of  the  optical  absorption  '  at  early  times  when  car¬ 
riers  are  assumed  to  be  in  a  thermal  equilibrium  with 
phonons,  i.e.,  TC=TL.  This  equation  can  also  be  regarded  as  a 
generalization  of  the  ABC  rate  model13  in  the  absence  of 
Shockley-Reed-Hall  nonradiative  recombination16  and  Auger 
recombination  for  defect-free  semiconductors  under  low  70. 
For  the  temporal  dependence,  we  assume  that  the  incident 
field  is  turned  on  at  time  t  =  0.  As  a  result,  the  initial  condi¬ 
tion  for  Eq.  (1)  is  simply  Nc( 0)  =  0,  where  we  have  neglected 
the  small  intrinsic  carrier  density  given  by  rij=spcpv 
Xexp(-EGl2kBTc)  with  pc=2(m*ekBTc/2irfi2)312  and  pv 
=  2{m*hkBTcl2Trh1)2'12  due  to  Nc>ni  at  low  temperatures. 
Here,  m*e  and  m*h  are  the  effective  masses  of  electrons  and 
holes,  respectively,  and  Ec  is  the  bare  band-gap  energy  of  the 
host  semiconductor  in  the  absence  of  photoexcited  carriers. 
We  will  look  for  the  steady- state  solution  of  Eq.  (1)  in  this 
paper. 

1.  Absorption  coefficient 

The  frequency-dependent  absorption  coefficient  in  Eq.  (1) 
is  given  by17 

/W1 «)  =  /?aL( w)  +  /3ab '  h  Im[aL(w)] ,  (2) 

nr(co)c 

where  eb  is  the  dielectric  constant  of  the  host  semiconductor 
and  aL(w)  =  af(to)  +  aeB(o))  is  the  total  Lorentz  function17  in¬ 
cluding  contributions  from  both  the  e-h  plasmas  (denoted  by 
the  superscript  pi),  as  well  as  from  the  excitons  (denoted  by 
the  superscript  ex).  Here,  the  scaled  refractive-index  function 
in  Eq.  (2)  is 

nr  M  =  ^={1  +  Re[aL(w)] 

+  V{1  +  Re[aL(W)]}2  +  {Im[ai(W)]}2}1/2.  (3) 
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The  Lorentz  function  for  the  photoexcited  e-h  plasmas  in 
the  long-wavelength  limit  is  calculated  as3,18 
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^pi 


2e2 


)(co  +  iy0 )  7 T( 


:0eb  ■>  0 


f 

Jo 


dk  k2\rvc(k)\2[l  —fe(k) 


-fhm 


Ee{k)  +  Eh(k) 


[h(w  +  iy0)f  -  \ Ee(k )  +  Eh{k)]2  ’ 


(4) 


where  hyQ  is  the  homogeneous  level  broadening  due  to  the 
finite  lifetime  of  quasiparticles,  k  is  the  wave  number  of  pho¬ 
toexcited  carriers,  and  the  plasma  frequency  of  the  e-h  com¬ 
posite  plasma  in  the  first  term  of  Eq.  (4)  is  defined  as 
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where  pte  and  nh  are  the  chemical  potentials  of  electrons  and 
holes.  For  given  Nc  and  Tc  at  each  moment,  the  chemical 
potential  /q  of  carriers  with  i  =  e  and  h  is  the  root  of  the 
following  equation: 
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In  Eq.  (4),  we  show  the  Lorentz  function  for  the  photoex¬ 
cited  e-h  plasmas.  The  Lorentz  function  for  the  .y-type  bound 
excitons  is  calculated  as4 
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where  fi,  is  the  reduced  e-h  mass.  In  addition,  the  interband 
dipole-moment  matrix  element  in  the  second  term  of  Eq.  (4) 
is 


rVc(£)|2  = 


E'g  h2 
4/r,.  (Ec  +  fi2k2l2pL,)2 


(6) 


The  kinetic  energy  Ej{k )  of  quasiparticles  for  i  =  e  or  h,  ini¬ 
tially  introduced  in  Eq.  (4),  in  the  screened  Hartree-Fock 
approximation  is  calculated  as4 


where  m0  is  the  free-electron  mass,  y'{)  is  the  dephasing  rate 
of  excitons,  n  =  1 , 2 ,  •  •  ■  is  the  radial  quantum  number  of 
i-type  bound  excitons,  En< 0  is  the  binding  energy,  and  aB 
=47 Te0ebh2l /jLre2  is  the  effective  Bohr  radius.  The  interband- 
transition  matrix  element  |PVC|2  in  Eq.  (11)  for  excitons  is 
given  by 
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Moreover,  we  define  rjn  in  Eq.  (11)  as  the  ratio  given  by 
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where  \k-k'\  =  x]k2  +  k’2-2kk'  cos  9,  0  is  the  angle  between  k 
and  k' ,  and  the  renormalized  energy  gap  in  Eq.  (6)  is  given 
by  E'c=Ee{0)+Eh(0).  The  last  term  in  Eq.  (7)  includes  the 
change  in  the  exchange  energy  due  to  screening  co2(q) 
=  Op[(l  +q2l q2)  +  {fiq2l 4-/j,r)2  and  the  static  screening  length 
1  lqs  in  the  Thomas-Fermi  limit  is  given  by17 
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where  the  exciton  eigenfunctions  (/'“(r)  and  eigenenergies  E„ 
in  Eq.  (11)  are  the  solutions  of  the  following  Schrodinger 
equation  for  a  screened  exciton: 


1  d 
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\  2/A 

r2  dr 

r~ 

l  dr  J 

4  t Te0ebr 


exp(-  qsr) 


=  0. 


(14) 


0“(r)  in  Eq.  (13)  denotes  the  bare  exciton  eigenfunction 
with  the  binding  energy  - EBln 1  where  EB=h2l2/j,raB.  In  ad¬ 
dition,  <fknx(r)  is  also  the  solution  of  the  above  Schrodinger 
equation  with  qs  =  0.  The  solution  of  Eq.  (14)  includes  the 
many-body  screening  effect,  i.e.,  E„  and  rjn  depend  on  qs  and 
decrease  with  increasing  carrier  density  or  carrier  tempera¬ 
ture.  A  more  extensive  treatment  of  nonlocal  excitons  has 
been  proposed  by  Rupper  et  al. 12 


The  result  in  Eq.  (7)  includes  the  many-body  energy  renor¬ 
malization  effect,  i.e.,  the  reduction  of  the  band-gap  energy 
Ec  with  increasing  carrier  density  or  with  decreasing  carrier 
temperature.  The  distributions  of  photoexcited  carriers,  ini¬ 
tially  introduced  in  Eq.  (4),  are  given  by  the  following  Fermi 
functions  for  i=e  and  h: 


m  = 


Ej{k )  -  fij 

kBTc  . 


(9) 


2.  Spontaneous  emission 

The  strong  interaction  of  e-h  plasmas  with  both  incident 
and  emitted  photons  forces  the  photon  and  charged  carrier 
systems  to  stay  in  a  thermal-equilibrium  state,  i.e.,  the  same 
temperature  for  both  photons  and  plasmas.  In  addition,  we 
know  that  the  radiative  lifetime  of  photoexcited  carriers  in 
the  host  semiconductor  is  known  to  be  on  the  order  of  1-10 
ns.  If  we  work  on  a  time  scale  slower  than  the  radiative 
lifetime  such  as  in  a  steady  state,  the  statistics  of  the  emitted 
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photons  will  obey  the  Bose-Einstein  distribution  at  each  mo¬ 
ment  with  the  same  temperature  as  that  of  the  carriers.3 

Based  on  the  assumption  that  charged  carriers  and  pho¬ 
tons  have  equal  temperature  Tc,  the  total  rate  of  spontaneous 
emission  per  unit  volume,  7£sp=7\!.j’p+7?.“  in  Eq.  (1),  can  be 
related  to  the  absorption  coefficient  /3abs[w)  through  the  so- 
called  Kubo-Martin-Schwinger  (KMS)  relation.19-22  Our  cal¬ 
culation  includes  the  level  broadening,  which  generalizes  the 
KMS  relation,  and  leads  to23 
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where  0(x)  is  the  step  function,  Xh  =  /xrl  m*h,  \e=firlm*e,  and 
\e+\h  =  I .  It  is  obvious  from  Eqs.  (15)  and  (16)  that  7*lsp 
will  be  enhanced  if  Nc  increases  since  it  is  proportional  to 
fJh-  Pp  could  also  be  enhanced  if  Ec  becomes  large  since  it 
is  proportional  to  <w3  or  E3C. 

B.  Dynamical  energy  balance  equation 

It  is  known  that  the  intraband  carrier-carrier  scattering 
time  is  on  the  order  of  100  fs,  which  is  much  faster  than  the 
radiative  lifetime  of  charged  carriers.  Therefore,  the  thermal- 
ization  of  hot  carriers  can  be  approximately  regarded  as  an 
adiabatic  process  compared  to  the  slow  radiative-decay  pro¬ 
cess,  i.e.,  the  total  energy  of  charged  particles  remains  con¬ 
served  at  each  moment.  Based  on  the  conservation  of  total 
energy  of  charged  particles,  the  dynamical  energy  equation 
for  the  e-h  plasmas  can  be  written  as7,9 

C0"~  =  /?abs(fiL)/0-  Wsp(Tc)  +  Wph(7077),  (17) 

where  the  heat  capacity  of  e-h  plasmas  is  given  by 
1  f00 

Q,  =  ,  2  dk  k^MEiik)  -  iitftm 

Z7TKb1  ci=e  h  J  o 

(IB) 

We  have  introduced  in  Eq.  (17)  the  power-loss  density  VV’ 

=  W9'  +  due  to  the  spontaneous  emission  of  photons.  yVpp 
and  y\^p  in  Eq.  (17)  are  calculated  by  inserting  the  factor  ha> 
into  Eqs.  (15)  and  (16),  respectively. 

When  both  the  defect  and  the  Auger  recombinations  are 
neglected,  the  steady-state  power  condition  can  only  be  sat¬ 
isfied  by  including  the  thermal  exchange  between  carriers 
and  phonons  as  shown  in  Eq.  (17).  However,  this  important 
thermal  exchange  term  has  been  dropped  in  the  papers  by 
Rupper  et  al.12  and  by  Sheik-Bahae  et  a/.13  when  they  used 
equal  carrier  and  lattice  temperature. 


The  initial  condition  for  Eq.  (17)  is  simply  Tc(0)  =  TL.  In 
this  paper,  we  will  look  for  the  steady-state  solution  of  Eq. 
(17).  Compared  to  the  previous  theory7,9  at  earlier  times,  we 
have  employed  the  generalized  KMS  relation  in  the  calcula¬ 
tions  of  WPp  and  VV('p .  Because  of  the  different  lattice  tem¬ 
perature  Tl  and  carrier  temperature  7’c,  there  exists  a  thermal 
energy  exchange  between  phonons  and  charged  carriers  in 
the  system,  which  is  represented  by  the  last  term  on  the 
right-hand  side  of  Eq.  (17)  and  is  calculated  by  using  the 
energy-dissipation  theorem  as9 
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where  \k+q\  =  (k2+q2  +  2kq  cos  0)112  and  the  photon-assisted 
phonon-scattering  process7  is  neglected  for  low  70.  In  Eq. 
(19),  we  assume  the  Bose  function  «ph(x)  =  l/[exp(x)-l]  for 
the  phonon  distribution,  fitoqK  is  the  phonon  energy  for  the 
wave  number  q  and  mode  X,  and  |C9\|2  represents  the 
carrier-phonon  coupling  matrix  element.  For  polar  semicon- 
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ductors  such  as  GaAs,  there  exist  both  acoustic  and  optical 
phonon  modes.9  For  optical  phonon  modes,  only  the 
longitudinal-optical  phonon  mode  strongly  couples  to  the 
charged  carriers.  Here,  we  set  \=LO,  LA,  and  TA  represent¬ 
ing  longitudinal-optical,  longitudinal-acoustic,  and  trans- 
verse-acoustic  phonons,  respectively.  In  addition,  we  get 
ojqk  =  <\q  for  \=LA  and  TA  by  applying  the  Debye  model24 
for  low-energy  acoustic  phonons  and  set  x=hcxq  in  the  sec¬ 
ond  term  of  Eq.  (19).  Using  the  Frohlich  model,24  we  find  the 
coupling  matrix  element  between  the  charged  carriers  and 
the  LO  phonons  in  Eq.  (19) 


\C, 


g,LOl 


^LoV  J_  _  1  j  e2 
2  /Woo  ej  e0(q2  +  q2)' 
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where  hu>L 0  is  the  energy  of  LO  phonons  and  ex  ( ev)  is  the 
high-frequency  (static)  dielectric  constant  of  the  host  semi¬ 
conductor.  For  the  acoustic -phonon  scattering  in  Eq.  (19),  on 
the  other  hand,  we  use  the  deformation-potential  approxi¬ 
mation.24  This  yields 
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where  cXA  and  cXA  are  the  sound  velocities  of  LA  and  TA 
phonons,  p0  is  the  ion  mass  density  of  the  host  semiconduc¬ 
tor,  D  is  the  deformation-potential  coefficient,  and  hu  is  the 
piezoelectric  constant.  In  our  calculations,  we  rewrite  Eq. 
(21)  using  q=x/hcLA  for  X  =  LA  and  Eq.  (22)  using  q 
=xlhc-iA  for  \  =  TA  corresponding  to  the  notation  in  the  sec¬ 
ond  term  of  Eq.  (19).  In  addition,  we  have  introduced  in  Eq. 
(19)  the  Debye  frequency  <xi^  =  (67r2A/a)1/3cx  for  acoustic 
phonons  and  the  maximum  momentum  transfer  fiqm 
=  \/4 iirtia>LO  for  optical  phonons,  where  Na  is  the  atom  num¬ 
ber  density.  The  screening  effect  on  the  carrier-phonon  scat¬ 
tering  is  included  in  Eqs.  (20)-(22)  in  the  Thomas-Fermi 
limit.  When  Tc  >  7),  we  find  VVph(7’r,  7)  )  <  0  from  Eq.  (19), 
indicating  the  loss  of  thermal  energy  from  hot  charged  carri¬ 
ers  to  cool  phonons.  For  a  fixed  lattice  temperature  TL,  the 
carrier  temperature  Tc(t)  can  be  found  from  Eq.  (17)  at  each 
time  t.  In  this  situation,  optical  carrier  cooling  (heating)  us¬ 
ing  nonthermal  e-h  plasmas  corresponds  to  tc<tl(tc 
>  Tf)  for  fixed  TL  (Ref.  9). 

Using  the  relation  between  Wsp  and  lZsp  and  the  steady- 
state  condition  of  Eq.  (1),  we  can  combine  the  first  two  terms 
on  the  right-hand  side  of  Eq.  (17)  into  one  term 
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FIG.  1.  (Color  online)  Schematic  for  the  shift  of  exciton  emis¬ 
sion  peak  [solid  (blue)  curves  in  the  upper  panel]  in  the  low-density 
regime  and  the  shift  of  e-h  plasma  emission  peak  [solid  (blue) 
curves  at  the  conduction-band  edge  denoted  by  rectangular  (cyan) 
blocks  in  the  lower  panel]  in  the  high-density  regime  with  increas¬ 
ing  laser  intensity  I0  [horizontal  (green)  dashed  arrows].  In  the  up¬ 
per  panel,  the  binding  energy  of  unscreened  excitons  is  -EBlrr  with 
n=  1,2,- while  the  binding  energy  of  screened  excitons  is  de¬ 
noted  by  En.  In  addition,  MlL  [upward  (brown)  arrow]  is  the  fixed 
laser  photon  energy.  In  the  lower  panel,  the  unrenormalized 
conduction-band  edge  is  Ec,  while  the  renormalized  conduction- 
band  edge  is  denoted  by  E'c.  As  explained  in  the  text,  a  very  weak 
switching  from  carrier  heating  to  carrier  cooling  is  expected  by 
increasing  I0  when  Nc  is  low  in  the  upper  panel,  while  an  opposite 
strong  switching  is  expected  when  Nc  is  high  in  the  lower  panel. 

In  general,  the  emission  spectrum  will  contain  two  peaks 
located  at  hu >=E'C+En  {En  <  0)  for  excitons  and  hu>=E'G  for 
e-h  plasmas.  Therefore,  if  the  laser  photon  energy  hil,  is 
below  the  emission-peak  energy,  we  expect  to  see  an  optical 
carrier  cooling  in  the  system.  If  the  laser  photon  energy  is 
above  the  emission-peak  energy,  on  the  other  hand,  we  will 
get  an  optical  carrier  heating.  As  a  result,  by  moving  hilL 
downward  across  E'c+En,  we  expect  to  see  the  crossover 
from  optical  carrier  heating  to  optical  carrier  cooling,  as  il¬ 
lustrated  in  the  upper  panel  of  Fig.  1,  if  the  exciton  emission 
dominates  when  I0  is  low  and  increases.  This  is  because  the 
exciton  binding  energy  En  decreases  with  70  or  Nc  in  nonsat- 
urated  absorption  regime  due  to  a  many-body  screening  ef¬ 
fect.  In  addition,  by  moving  /HI,  upward  across  E'a,  we  will 
find  the  change  from  optical  carrier  cooling  to  optical  carrier 
heating,  as  shown  in  the  lower  panel  of  Fig.  1,  if  the  plasma 
emission  dominates  when  /0  is  high  and  increases.  This  is 
because  the  band-gap  reduction  EG-E'G  increases  with  I0  due 
to  the  many-body  energy  renormalization  effect.  In  reality, 
however,  the  optical  carrier  cooling  at  low  70  is  too  weak  to 
be  seen  since  the  emission  peak  in  dlZsp(w)l dco  becomes 
negligible  at  very  low  Nc. 


TIT.  NUMERICAL  RESULTS  AND  DISCUSSIONS 


/Sabs(D/J/0-Wsp(7’c)  =  J  dw 

0  ait  liiLtu  as  iunuw  s.  ni(),  lit;,  —  vi.uz.  iiiq,  i^G 

(23)  =  1.519-5.405  X10-4[7’2/(7’L+204)](eV-K-1),  ev=  13.18, 

600=10.89,  eb=( 65  +  600)/ 2,  hcjLO=  36.25  meV,  /rwx0 
=  33.29  meV,  cLA=5.14X  105  cm/s,  cxa=3.04X  105  cm/s, 
Po  =  5.3  g/cm3,  D  =  -9.3  eV,  /?14=1.2X107  V/cm,  Na 


dK^w) 

dco 


In  our  numerical  calculations,  we  have  chosen  GaAs  as 
(hilL  -  hw) .  the  host  semiconductor.  The  parameters  for  this  host  material 


Here,  the  positive  d'IZsp( oj)  /  doj  represents  the  emission  spec¬ 
trum  as  a  function  of  the  energy  hu>  of  emitted  photons. 
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ha-  Eq  (  meV  ) 
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FIG.  2.  (Color  online)  Absorption  spectrum  /8abs(w)  [in  (a)  with 
a  unit  of  104  cm-1]  and  spontaneous-emission  spectrum 
dlZsp(co)l  da  [in  (b)  with  an  arbitrary  unit  (a.u.)]  as  functions  of  the 
photon  energy  ha  with  respect  to  EG.  Here,  the  dashed  (blue) 
curves  correspond  to  a  lower  laser  intensity  70=10-3  W/cm2  ( Nc 
=  3.41  X  1010  cm-3),  while  the  solid  (red)  curves  are  associated 
with  a  higher  laser  intensity  /0=103  W/cm2  (Nc=9.55 
X  1015  cm-3).  The  label  “X0.001”  in  (b)  indicates  that  the  peak  of 
the  dashed  (blue)  curve  has  been  increased  by  103  times. 

=4.42  X  1022  cm-3,  h y0= h yQ= 0. 05 EB= 0.21  meV,  and  TL 
=4  K.  The  above  expression  for  Ec  reflects  the  fact  of  the 
reduction  of  the  bare  band-gap  energy  with  increasing  lattice 
temperature.  The  other  adjustable  parameters  such  as  70  and 
hHL  will  be  directly  given  on  the  figures. 

We  present  in  Fig.  2  the  calculated  steady-state  absorption 
coefficient  ySabs(w)  from  Eq.  (2)  [in  (a)],  as  well  as  the  emis¬ 
sion  spectrum  dRsp(a)/ da  from  Eqs.  (15)  and  (16)  [in  (b)]  as 
a  function  of  the  probe  photon  energy  ha-Ec  at  fiLlL- Ec 
=  -3.8  meV  for  70=10-3  W/cm2  (blue  curves)  and  70 
=  103  W/cm2  (red  curves),  respectively.  From  Fig.  2(a)  at 
70=  10-3  W/cm2,  we  see  a  very  strong  absorption  peak  from 
the  n=  1  exciton  state  below  the  conduction-band  edge,  as 
predicted  by  Eq.  (11),  in  addition  to  another  two  weak  ab¬ 
sorption  peaks  from  n= 2  and  n  =  3  exciton  states.  However, 
the  step  feature  in  /3abs(<w)  due  to  absorption  by  e-h  plasmas, 
as  predicted  by  Eq.  (4),  at  the  conduction-band  edge  is  al¬ 
most  unresolved  in  this  case.  When  the  incident  laser  inten¬ 


FIG.  3.  (Color  online)  Ratio  of  the  calculated  carrier  tempera¬ 
ture  Tc  to  the  fixed  lattice  temperature  TL  [squares  (blue)  curve  and 
in  a  left  scale]  and  the  calculated  density  of  photocarriers  [triangles 
(red)  curve  and  in  a  right  scale  with  a  unit  of  10 15  cm-3]  as  func¬ 
tions  of  the  laser  intensity  70  (with  a  unit  of  kW/cm2)  at  TL=4  K 
for  MlL-EG=- 3.8  meV. 

sity  70  is  high,  a  large  density  of  photoexcited  carriers  effec¬ 
tively  screens  the  exciton  attractive  interaction,  as  predicted 
in  Eq.  (14)  with  a  large  value  for  qs.  As  a  result,  the  exciton 
binding  energy  becomes  negligibly  small.  At  the  same  time, 
the  ratio  of  the  dipole  moments  for  interband  transitions  in 
Eq.  (13)  is  also  greatly  suppressed.  These  two  factors  com¬ 
bined  together  are  responsible  for  the  disappearance  of  the 
exciton  absorption  peaks  at  70=103  W/cm2.  In  this  case, 
however,  the  step  feature  in  /3abs(«)  at  the  conduction-band 
edge  is  enhanced  due  to  the  suppressed  background  from  the 
exciton  absorption  peaks.  The  exciton  and  plasma  effects  can 
also  be  seen  from  the  emission  spectrum  in  Fig.  2(b).  For 
70=10-3  W/cm2,  a  weak  emission  peak  occurs  below  the 
conduction-band  edge  for  the  n  =  1  exciton  state,  which  is 
replaced  by  a  strong  peak  from  the  e-h  plasmas  at  the 
conduction-band  edge  when  70  =  1 03  W/cm2.  Moreover,  the 
peak  value  in  dRsp(a)l  da  is  found  to  decrease  with  Tc. 

Figure  3  displays  the  calculated  steady-state  carrier  tem¬ 
perature  TcITl  (blue  curve)  and  density  Nc  of  photoexcited 
carriers  (red  curve)  as  functions  of  the  incident  laser  intensity 
Iq  for  hilL-Ec=-3.8  meV.  In  this  case,  the  laser  pumping 
is  set  under  the  conduction-band  edge  since  h£lL<EG,  and 
the  exciton  effect  is  negligible  for  the  values  of  70  in  the 
shown  range.  When  70  is  increased  above  2  kW/cm2,  we 
find  that  Nc  determined  from  Eq.  (1)  becomes  nearly  inde¬ 
pendent  of  Iq  due  to  saturated  optical  absorption.3  At  the 
same  time,  however,  TCITL  determined  from  Eq.  (17)  still 
drops  as  70  increases,  indicating  optical  carrier  cooling  due  to 
TcITl<1  in  contrast  to  electrical  carrier  cooling.2  The  in¬ 
creasing  Nc  with  I0  will  shrink  the  band  gap  from  Ec  to  E'c 
determined  by  Eq.  (7),  where  the  renormalized  conduction- 
band  edge  E'c  roughly  labels  the  peak  energy  in  emission 
spectrum  dRsJa)l  da.  It  is  evident  that  TCITL  will  eventually 
increase  with  70  due  to  optical  carrier  heating  after  passing 
through  a  minimum,  as  predicted  by  Eqs.  (17)  and  (23),  simi¬ 
lar  to  electrical  carrier  heating.2 

We  show  in  Fig.  4  the  calculated  steady-state  carrier  tem¬ 
perature  TcITl  (blue  curve)  and  peak  energy  fiapX-EG  of 
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FIG.  4.  (Color  online)  Ratio  of  the  calculated  carrier  tempera¬ 
ture  Tc  to  the  fixed  lattice  temperature  TL  [squares  (blue)  curve  and 
in  a  left  scale]  and  the  peak  energy  hwp\-EG  of  emission  spectrum 
with  respect  to  the  bare  band-gap  energy  Ec  [triangles  (red)  curve 
and  in  a  right  scale  with  a  unit  of  meV]  as  functions  of  the  laser 
intensity  70  (with  a  unit  of  W  /  cm2  and  a  logarithmic  scale)  at  TL 
=4  K  for  mL-EG=-3.S  meV. 
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FIG.  5.  (Color  online)  Difference  of  the  calculated  carrier  tem¬ 
perature  Tc  and  the  fixed  lattice  temperature  TL  [squares  (blue) 
curve  and  in  a  left  scale  with  a  unit  of  K]  and  the  calculated  density 
of  photocarriers  [triangles  (red)  curve  and  in  a  right  scale  with  a 
unit  of  1015  cm-3]  as  functions  of  the  detuning  of  laser  photon 
energy  hClL-Ec  with  respect  to  the  bare  band-gap  energy  Ec  (with 
a  unit  of  meV)  at  77=4  K  for  70  =  5  kW/cm2. 


photoluminescence  (red  curve)  as  functions  of  the  incident 
laser  intensity  70  for  hilL—EG=5  meV  (the  laser  pumping 
energy  is  set  above  the  conduction-band  edge  due  to  /HI, 
-Eg> 0).  When  70  is  increased  up  to  1  W/cm2,  the  exciton 
effect  dominates  in  this  range  of  70.  Therefore,  hu>pi-Ec 
shifts  up  from  -4.4  to  -3.2  meV  due  to  the  screening  effect 
on  the  exciton  binding  energy.  When  70  is  increased  from 
10  W/cm2,  the  emission  process  will  be  dominated  by  e-h 
plasmas.  In  this  case,  ha>pl-Ec  scales  with  the  Fermi  edge  of 
e-h  plasmas,  which  moves  up  in  energy  with  Nc  or  equiva¬ 
lently  with  Iq.  When  70  changes  between  1  and  10  W/cm2, 
the  emission  process  switches  from  exciton  dominance  to 
e-h  plasma  dominance  in  this  crossover  region.  Correspond¬ 
ingly,  TcITl  remains  unity  before  the  e-h  plasma  dominance 
begins.  However,  TCITL  increases  dramatically  with  70  once 
the  emission  process  is  dominated  by  e-h  plasmas,  indicating 
a  strong  optical  carrier  heating  in  the  system. 

Figure  5  exhibits  the  calculated  steady-state  temperature 
difference  TC—TL  (blue  curve)  and  density  Nc  of  photoex- 
cited  carriers  (red  curve)  as  functions  of  the  frequency  de¬ 
tuning  hClL-Ec  of  a  pump  laser  at  70=5  kW/cm2.  For  such 
a  large  value  of  70,  the  e-h  plasmas  fully  dominate  the  optical 
process.  In  the  range  of  hElL>EG,  Nc  increases  with  the 
frequency  detuning  in  a  superlinear  dependence  since  the 
Fermi  energy  is  proportional  to  the  frequency  detuning  due 
to  suppressed  saturation  for  optical  absorption.  On  the  other 
hand,  in  the  opposite  range  of  hD,L<EG,  Nc  decreases  with 
the  frequency  detuning  due  to  the  enhanced  peak  value  in 
dRsp(w)l da>,  and  it  reaches  a  minimum  slightly  below  the 
band  gap  Ec.  Moreover,  Tc  increases  with  the  frequency  de¬ 
tuning  from  below  TL  in  the  range  of  /HI,  -  eg<  0  when 
hHL  approaches  the  emission-peak  energy  around  the 
conduction-band  edge,  as  can  be  seen  from  Eq.  (23).  Once 
/H  I,  >  Ec  is  reached,  Tc  is  initially  locked  to  7’,  and  becomes 
independent  of  the  frequency  detuning  due  to  saturated  opti¬ 


cal  absorption.  This  is  followed  by  a  very  strong  increase  in 
Tc  with  the  frequency  detuning  due  to  enhanced  optical  ab¬ 
sorption  after  the  saturation  is  suppressed.  Finally,  the  in¬ 
crease  in  Tc  is  greatly  slowed  down  due  to  a  reduced  peak 
value  in  dRsp{(o)  /  d(o  with  increasing  Tc.  This  demonstrates  a 
switching  from  optical  carrier  cooling  to  heating  when  the 
frequency  detuning  of  a  pump  laser  sweeps  from  below  the 
band  gap  to  above  the  band  gap. 

IV.  CONCLUSION 

In  conclusion,  we  have  applied  the  dynamical  density  and 
energy  equations  to  study  the  steady-state  optical  carrier 
cooling  in  intrinsic  semiconductors  mediated  by  nonthermal 
electron-hole  plasmas,  where  the  lattice  system  is  assumed  in 
thermal  contact  with  an  external  heat  bath.  We  have  included 
many-body  effects  such  as  Coulomb  renormalization  of 
band-gap  energy,  screening,  and  excitonic  interaction  in  the 
calculations  of  optical  absorption,  spontaneous-emission  and 
carrier-phonon  scattering.  We  have  demonstrated  and  ex¬ 
plained  the  many-body  effects  on  optical  carrier  cooling  in 
the  system  at  low  lattice  temperatures  and  studied  the  depen¬ 
dence  of  optical  carrier  cooling  on  the  pump-laser  intensity 
and  photon  energy. 

Auger  recombination  has  not  been  included  in  our  calcu¬ 
lations  for  low  lattice  temperatures  and  low  photocarrier  den¬ 
sities.  We  would  like  to  point  out  that  in  the  current  case,  the 
weak  Auger  recombination  in  our  system  is  expected  to 
cause  a  slight  modification  to  the  steady-state  photocarrier 
density  determined  from  Eq.  (1)  and  a  heating  of  the  lattice 
through  a  multiphonon  emission  process.  The  modification 
of  the  photocarrier  density  will  not  alter  the  predicted  many- 
body  effects  on  optical  carrier  cooling  qualitatively.  In  addi¬ 
tion,  the  heat  dissipated  to  the  lattice  by  a  left-over  high- 
energy  excited-state  carrier  in  a  weak  Auger  recombination 
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process  will  be  completely  absorbed  by  the  external  heat 
bath  to  maintain  a  fixed  lattice  temperature. 

With  a  similar  procedure  used  by  Rupper  et  al.}2  we  will 
calculate  the  following  ratio  to  characterize  the  importance 
of  the  Auger  recombination  with  respect  to  the  radiative  re¬ 
combination,  which  is  given  by 


Auger-recombination  coefficient  and  B(T)  is  the  radiative- 
recombination  coefficient  at  the  temperature  T.  Using  T 
=  10  K,  Nc=  1016  cm“3,  and  5(10  K)  =  3X10"4  cm3/s,  we 
find  [C(10  K)NC/B(10  K)]<S1,  which  justifies  neglecting 
the  Auger  recombination  in  this  paper. 


C{T)N 3  C(300  K) 

B{T)N]  ~  B(T)  £XP 


2.24  1  - 


300  \ 
T(K)) 


Nc{  cm"3), 


where  C(300  K)=4X  10  30  cm6/s  is  the  room-temperature 
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